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We have calculated the splitting functions governing the evolution of the un- 
polarized parton distributions of the photon at the next-to-next-to-leading 
order (NNLO) of massless perturbative QCD. The results, presented here 
mainly in terms of compact and accurate parametrizations, are consis- 
tent with our previous approximations based on the lowest six even-integer 
Mellin moments. Consequently the NNLO corrections are small in both the 
MS and the DIS-y factorization schemes at momentum fractions x > 0.1. 



PACS numbers: 12.38.Bx; 13.60.Hb, 13.66.Bc, 14.70.Bh 



The partonic structure of the photon, accessed in particular by the deep- 
inelastic structure function F2 , is a classic subject in perturbative QCD. 
The leading-order (LO) and next-to-leading order (NLO) expressions for 
the photon-parton splitting functions Pp-y, p = q, g, and the coefficient 
functions C2^^ have been known for a long time [1-6]. A couple of years 
ago we have presented [7] the corresponding next-to-next-to-Ieading order 
(NNLO) corrections, albeit with one important qualification concerning the 
©(aOg) splitting functions: at that time we were only able to derive the 
lowest six even- integer Mellin moments, = 2, . . . , 12, which are sufficient 
for a reliable reconstruction only at momentum fractions x 0.05. 
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Using methods and results derived in the meantime for the hadronic 
case [8-10], we are now able to present the complete results for the splitting 
functions and coefficient functions to order aOg, thus finalizing the NNLO 
description and providing an important partial result for F2 at N^LO. In 
the present brief contribution, we confine ourselves to the NNLO splitting 
functions, presenting accurate compact parametrizations of those results 
which are rather lengthy. Furthermore we take a first look at their numerical 
effects. A more detailed account, including the exact results and the aa^ 
photonic coefficient functions for F2 and Fl, will be presented elsewhere [11]. 

At lowest order in the electromagnetic coupling = aem/(4vr), the 
parton distributions of the photon are subject to the evolution equations of 
the form 

-^ = P^ + P0Q^, (1) 

where ^ represents the MS factorization and renormalization scale, and ® 
stands for the Mellin convolution in the momentum variable. For brevity 
writing out only the flavour-singlet case, is given by 
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n£ is the number of active flavours, and the splitting-function matrices are 

= {^^A . P = ( ) . (3) 

The expansions of the photon-parton and parton-parton splitting functions 
up to NNLO read, with Og = as/(4vr), 

P^ = (pW+a^pW+a^pf) (4) 

P = asP(°) +alP^^'^ + alP^^^ . (5) 

The hadronic quantities P^^^ in Eq. © can be found in Refs. [8,9]. Our 
new results for the photonic non-singlet (see Ref. [7] for notational details) 
and gluon MS splitting functions P^^^ in Eq. Q can be parametrized as 

KsP^ii^) - 128/27 112/9 175.3 + 142.3 Li + 1353 

- 1262 X + 449.2x2 - 1445 - LqLi (162.7 Lq + 195.4 Li) + 1169 xLq 
+ 50.08 (1 - x)Ll + 744.6 Lq + 201.6 + 80/3 + 64/27 

+ nf{ - 32/27 Lf - 11.858 - 18.77 Li - 40.035 + 114.4 x 

- 24.86x2 - 53.39x3 _^ ^^^^ (8.523 Lq + 269.4 Li) - 26.63 xLq 

+ 270.0 (1 - x)L\ - 21.55 Lq - 10.992 - 32/27 } (6) 
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and 

5-^P^^\x) ^ (1 - x) { 32/27 - 79.13 L? + 87.22 Li + 1738 - 1580 X 

- 160.0x2 - 566.7x^ - LqLi (549.5 + 1230 Lq + 433.2 Li) + 2176 Lq 
+ 1123.7 Ll + 2400/27 + 448/27 - 73.1409 x ~^ + 128/3 x ^^Lq} 

+ rif (1 - x) { - 32/9 Ll + 16.38 Li + 68.10 - 36.42 x + 56.95 x^ 

- 44.10 x^ - LqLi (16.18 + 38.33 Lq + 9.133 Li) - 10.76 Lq 

+ 26.41Lg- 64/27 Lg-40.5597x"^} (7) 

with Lq = Inx and Li = ln(l — x). These parametrizations deviate from 
the lengthy full expressions by about 0.1% or less, an accuracy which should 
be amply sufficient for practical applications. On the other hand, the exact 
expression is very compact for the NNLO pure-singlet splitting function. 



S-^P^'^lix) = 4/3 nf { 2464/81 x"^- 432 - 72 x + 38360/81 x2 

- Lo (344 + 368x + 3584/27 x^) - (144 + 104x + 224/9 x^) 

- Ll (16 - 16 X - 128/9 x^) - L^ 8/3 (1 - 2 x)} . 
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Recall that Pq^ in Eqs. (jSJ and @ is given by Pq^ = 5s/<^ns -fns7 + -fps7 with 
(^s = 3 J2jLi ■ A pure-singlet term enters at order aa^ for the first time. 
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Fig. 1. The exact results for the acmO-s photon-quark (left) and photon-gluon 
(right) splitting functions (multiplied by x) in the MS scheme, compared with the 
(dashed) estimated error bands based on the lowest six even- integer moments [7]. 
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Fig. 2. The perturbative expansion Q of the photon-quark (left) and photon- 
gluon (right, muhiphed by a;) sphtting functions in the DIS-y factorization scheme 
for typical values of as and rif. Note the very different scales of the two graphs. 

Our complete results for Pp-y , p = q, g, are compared in Fig. ^ with the 
previous approximations of Ref. [7]. In Fig. [51 these results are combined, 
after transformation to the DIS^ scheme [6,7], with the lower-order splitting 
functions. As indicated by the previous approximate results, the perturba- 
tive expansion is well-behaved at least at x > 0.1. For a further discussion, 
including the small-x limit, the reader is referred to Ref. [11]. 
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